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ABSTRACT 

Subject of this paper are the statistical properties of ellipticity alignments between galaxies 
evoked by their coupled angular momenta. Starting from physical angular momentum mod- 
els, we bridge the gap towards ellipticity correlations, ellipticity spectra and derived quantities 
such as aperture moments, comparing the intrinsic signals with those generated by gravi- 
tational lensing, with the projected galaxy sample of EUCLID in mind. We investigate the 
dependence of intrinsic ellipticity correlations on cosmological parameters and show that in- 
trinsic ellipticity correlations give rise to non-Gaussian likelihoods as a result of nonlinear 
functional dependencies. Comparing intrinsic ellipticity spectra to weak lensing spectra we 
quantify the magnitude of their contaminating effect on the estimation of cosmological param- 
eters and find that biases on dark energy parameters are very small in an angular-momentum 
based model in contrast to the linear alignment model commonly used. Finally, we quantify 
whether intrinsic ellipticities can be measured in the presence of the much stronger weak lens- 
ing induced ellipticity correlations, if prior knowledge on a cosmological model is assumed. 
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1 INTRODUCTION 

Weak cosmic shear, i.e. lensing by the gravitational field of the cos- 
mic matter distribution (Blandford et al. 1991; Seitz et al. 1994; 
Seitz & Schneider 1994; Kamionkowski et al. 1998), is considered 
to be an excellent probe of structure formation processes, preci- 
sion measurements of cosmological parameters (Hu 1999, 2002a,b; 
Takada & White 2004; Hannestad et al. 2006) and the influence 
of dark energy on cosmic structure formation (Huterer & Turner 
2001; Huterer 2002, 2010; Amara & Kitching 2011; Kunz 2012). 
The primary observable are ellipticity correlation functions or their 
Fourier counterparts (Jain & Seljak 1997; Hu & Tegmark 1999; Hu 
& White 2001; Hu & Jain 2004). These shape correlations have 
been first detected by a number of research groups more than 10 
years ago (Van Waerbeke et al. 2000; Kaiser et al. 2000; Bacon 
et al. 2000; Wittman et al. 2000) and are now routinely used for 
parameter estimation. Correlations in shapes of galaxies are intro- 
duced because light rays from neighbouring galaxies experience 
correlation distortions due to correlations in the tidal fields through 
which the respective rays propagate. A common assumption is the 
absence of intrinsic correlations such that any positive shape cor- 
relation can be attributed to the gravitational lensing effect. This 
hypothesis, however, might be flawed as there are physical mecha- 
nisms by which galaxies are intrinsically shape correlated: Due to 
the fact that neighbouring galaxies form from correlated initial con- 
ditions, their respective angular momenta are correlated (Croft & 
Metzler 2000; Heavens et al. 2000; Crittenden et al. 2001; Mackey 
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et al. 2002). Assuming that the galactic disks are established with 
their symmetry axes colinear with the host haloes' angular momen- 
tum directions one would observe galaxies at correlated angles of 
inclination and therefore with correlated ellipticities. 

A possible consequence of this new source of ellipticity cor- 
relation is its interference with the determination of cosmological 
parameters from weak lensing data, in particular the properties of 
dark energy. This issue has been the target of a number of investi- 
gations: Commonly, the description of intrinsic ellipticity correla- 
tions was based on the linear alignment model (Catelan et al. 2001; 
Hirata & Seljak 2004), 



C (d 2 x - <9j) O and e x = 2C d x d y <b. 



(1) 



which provides a direct modelling of the ellipticity field on the tidal 
shears d a dp(b (here, the z-axis of the coordinate system is aligned 
with the line-of-sight) and is able to give a consistent description of 
gradient and vorticity modes of the ellipticity field. The constant of 
proportionality was fixed by comparison with observations (Bridle 
& King 2007; Joachimi et al. 2012). 

If galaxy ellipticities are in fact described by an alignment 
model linear in the tidal fields, cosmological parameters, in partic- 
ular the dark energy equation of state parameters would be severly 
biased (Bridle & King 2007; Joachimi & Bridle 2010; Kirk et al. 
2010, 2011). Apart from ellipticity correlations themselves, ellip- 
ticity position-correlations were affected and ellipticity data would 
exhibit cross-correlations between intrinsic ellipticities and weak 
lensing (see, in particular, Hirata & Seljak 2004). 

There are basically four ways of dealing with intrinsic alig- 
ments. Firstly, they can be removed from data by using the fact that 
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they are a small scale phenomenon (Hey mans & Heavens 2003; 
King 2005; King & Schneider 2002, 2003) which takes place at 
the cost of increasing statistical uncertainties. Secondly, one can 
take advantage of the fact that intrinsic alignments have different 
statistical properties in comparison to weak lensing ellipticity cor- 
relations (Crittenden et al. 2002; King & Schneider 2003), most 
notably it is possible to use the statistics of vortical excitations in 
the ellipticity field which are exclusively sourced by intrinsic align- 
ments. Thirdly, one can design line-of-sight weightings that null 
out contributions due to intrinsic alignments (King 2005; Joachimi 
& Schneider 2009, 2008) which marginally increase statistical un- 
certainties on cosmological parameters. Parameter inference from 
spectra that result from data in this way still yields unbiased esti- 
mates. Finally, one can parameterise the intrinsic alignment contri- 
bution to weak lensing data and have those model parameters be 
determined by data alongside the cosmological parameters under 
consideration. Maginalisation over the parameters entering the in- 
trinsic alignment model then propagates the statistical errors of the 
alignment model on to the cosmological model. With a physically 
correct alignment model the estimates of cosmological parameters 
will remain unbiased. The feasibility of this approach under the as- 
sumption of Gaussian likelihoods has been demonstrated (Bridle & 
King 2007; Bernstein 2009; Joachimi & Bridle 2010; Kirk et al. 
2010; Laszlo et al. 2012). 

The motivation of this work was to explore intrinsic alignment 
effects and their observable properties in angular momentum-based 
alignment models. In these models, the ellipticity is quadratic in 
the tidal shear field and because they use in principle a mechani- 
cal model of angular momentum generation and ellipticity align- 
ment, the model parameters can be constrained from information 
other than ellipticity data. We will need two physically meaning- 
ful variables: a parameter which is related to the angular momen- 
tum model and whose value can be measured in structure forma- 
tion simulations and a disk morphology parameter which is acces- 
sible in galaxy surveys. Clearly, quadratic alignment models will 
differ in their prediction of ellipticity correlations compared to lin- 
ear alignment models. Together with the above results employing 
linear alignment models we hope to complete the view on intrinsic 
alignments and their relevance for future weak lensing surveys. 

The aim of this paper is threefold: (i) We investigate and com- 
pare two angular-momentum based alignment models in their pre- 
dictions for ellipticity correlations and formulate these predictions 
in terms of ellipticity correlation functions, ellipticity spectra and 
the scale-dependence of the variance of the ellipticity fields and 
compare these predictions with the equivalent quantities sourced 
by weak gravitational lensing (Sects. 2 and 3). (ii) The dependence 
of the two ellipticity models in consideration on cosmological pa- 
rameters is investigated and their likelihoods are derived. With this 
knowledge, we quantify the contamination of weak lensing data 
with an intrinsic alignment contribution and quantify how this con- 
taminations impacts on the estimation of cosmological parameters 
(Sect. 4). (iii) We investigate if there is a possibility of observ- 
ing intrinsic correlations in the presence of much stronger lensing- 
induced ellipticity correlations and develop statistical methods for 
answering these questions (Sect. 5). Throughout we will focus on 
intrinsic ellipticity correlations caused by correlated angular mo- 
menta, which is an applicable model for spiral galaxies. Those 
intrinsic alignments are proportional to the squared tidal field, in 
contrast to the linear alignment model valid for elliptical galaxies. 
In this limit, we neglect cross-correlations between intrinsic ellip- 
ticity alignments with the tidal field and gravitational lensing, as 
those correlations are proportional to the expectation value of the 



tidal field cubed, which vanishes in the case of Gaussian statistics. 
Specifically, we consider the case of EUCLID'S weak lensing sur- 
vey (Amendola et al. 2012). 

The reference cosmological model used is a spatially flat 
vi'CDM cosmology with Gaussian adiabatic initial perturbations 
in the cold dark matter density field. The parameter choice is 
motivated by the WMAP7 results (Komatsu et al. 2011; Larson 
et al. 2011): Sl m = 0.25, n s = 1, <x 8 = 0.8, Cl b = 0.04 and 
H = 100/jkm/s/Mpc, with h = 0.72. The dark energy equation of 
state is set tow = -0.95. 



2 COSMOLOGY 

2.1 Dark energy cosmologies 

In spatially flat Friedmann-Lemaitre Robertson- Walker cosmolo- 
gies with the matter density parameter Q„, and a dark energy com- 
ponent with equation of state w(a), the Hubble function H(a) = 
d In a I At is given by 

H 2 (a) a m I r' \ 

= + (1 - OJ exp ^3 J d In a (1 + w(a))j . (2) 

The value w = -1 corresponds to the cosmological constant A. The 
Hubble function describes the time evolution of the metric and can 
be used for relating comoving distance ;f and scale factor a: 

X = c f da ^— , (3) 

in units of the Hubble distance^ = c/H . The Hubble function 
also determines the critical density, p cr jt = 3H 2 /(&nG). 



2.2 CDM power spectrum 

The linear CDM density power spectrum P(k) describes the fluctu- 
ation amplitude of the Gaussian homogeneous density field 6, 

(6(k)5(k')) = (2nf6 D (k + k')P(k), (4) 

and is given by the ansatz P(k) oc k" s T 2 (k) with the transfer function 
T(k). In cosmologies with low Cl m , T(k) is fitted by (Bardeen et al. 
1986; Sugiyama 1995): 

T(q) = ln(1 2 + 34 g 34t?) + 3 - 89fl + t 16 - 1 ^ + < 5 - 46 <7) 3 + (6.71g) 4 ) 4 .(5) 

The wave number k is rescaled with the shape parameter T Si m h, 
q = k/T. The spectrum P{k) is normalised to the variance cr 8 of the 
density field on scales of R = 8 Mpc /h, 

< r l = ^f dkk 2 P(k)W 2 (kR) = J din it A 2 (k)W 2 (kR). (6) 

W(x) = 3ji(x)/x is the Fourier-transformed spherical top hat filter 
function. je(x) refers to the spherical Bessel function of the first 
kind of order I (Abramowitz & Stegun 1972; Arfken & Weber 
2005) and the dimensionless variance per logarithmic wavenum- 
ber A 2 (k) = k 3 P(k)/(27T 2 ) can be used instead of the CDM spec- 
trum P(k). In computing ellipticity correlation functions and ellip- 
ticity spectra we will employ a smoothed CDM spectrum P(k) — » 
P (k) exp(-(kR) 2 ) with a smooothing scale R that corresponds to a 
mass cutoff at a halo mass M. Those two quantities are related by 
M = 4tt/3 fi,„p crit R\ 
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2.3 Linear structure growth 

As long as the amplitudes in the cosmic density field are small, 
5 <k 1, the density field grows in a homogenous way, S(x, a) = 
D + (a)6(x,a = 1). The growth function D + (a) results from solving 
the growth equation (Turner & White 1997; Wang & Steinhardt 
1998; Linder & Jenkins 2003), 



1 



— D + (a)+- 3 + 
da 1 a\ a ma 



dmH\ d 



da 



D + (a) 



2a 1 



Q. m (a)D + (a). 



(7) 



Nonlinear structure formation enhances the CDM-spectrum P(k, a) 
on small scales by one and a half order of magnitude, which is 
described by the fit suggested by Smith et al. (2003). 



2.4 Angular momenta from tidal shearing 

Angular momenta of dark matter haloes embedded in potential 
flows in the large-scale structure are generated by a mechanism ref- 
ered to as tidal shearing, where the differential motion of mass ele- 
ments inside a protohalo gives rise to a torquing moment (Hoyle 
1949; Sciama 1955; Peebles 1969; Doroshkevich 1970; White 
1984): 



a s H{a) 



dD + 
da 



(8) 



i.e. it is the variation d a vp of the velocities vp ~ <9g© across the 
protohalo and hence the tidal field <t> a p 



4> 



d a d p <b 



(9) 



which is responsible for angular momentum generation. The mass 
distribution inside the protohalo itself is described by its inertia ten- 
sor /. 



aj3> 



up 



t a 3 f d 3 



q(q- qh(q - q)p 



(10) 



i.e. the second moments of the matter distribution, with the centre 
of mass at the position q and the integration comprising the La- 
grangian volume of the protohalo. Throughout, we use Einstein's 
summation convention. 

The Levi-Civita-symbol in eqn. (8) generates the interesting 
misalignment property between the shear and inertia eigensys- 
tems which is required for generating angular momentum: Only 
the antisymmetric contribution XjL, derived from the commutator 



fty 



[Ipy, 0/3 r ], to the product between the tensors /g r and is 
non- vanishing in contraction with the antisymmetric e a p y and there- 
fore relevant for angular momentum generation. The symmetric 
contribution Xt, which can be isolated using the anticommutator 



U/iy, <&py} cancels in the contraction. This means that for an- 



gular momentum build-up, the tidal shear and the inertia tensors are 
not allowed to have a common eigensystem and be skewed relative 
to each other (Schafer 2009; Schafer & Merkel 2012). Likewise, 
degeneracies in this relation evoked by spatial symmetries in the 
two tensors can prohibit the generation of angular momentum. Al- 
ternative models of galaxy angular momenta assume that the haloes 
are spun up by non-central, anisotropic infall in filaments (see, for 
instance, Pichon et al. 2010; Kimm et al. 201 1; Codis et al. 2012). 



2.5 Galaxy ellipticities 

Ellipticity correlations between galaxies are traced back to corre- 
lated angular momenta of their host haloes. CDM haloes acquire 
their angular momentum by tidal shearing and due to the fact that 



neighbouring galaxies experience correlated tidal fields, their an- 
gular momenta are correlated in consequence. The direction of the 
angular momentum L in turn determines the angle of inclination 
under which the galactic disk is viewed, and ultimately the elliptic- 
ity which is attributed to the galactic disk (Heavens et al. 2000; Crit- 
tenden et al. 2001, 2002; Mackey et al. 2002; Heymans & Heavens 
2003): Linking the angular momentum direction L = L/L to the 
components of the complex ellipticity e using the above argument 
yields 



with 



'- 7— . 

1 +L? 



2a 



1 +LV 



(ID 



if the coordinate system is aligned with its z-axis being parallel to 
the line of sight. A rotation of the coordinate frame by tp causes the 
complex ellipticity to rotate twice as fast, e — > exp(2i<^)e. a is a 
free parameter weakening the dependence between inclination an- 
gle and ellipticity for thick galactic disks and has been determined 
to be (x — 0.75 in the APM sample (Crittenden et al. 2001) with a 
large uncertainty. 

It should be emphasised that the assumption of a galactic disk 
forming perpendicularly to the host halo angular momentum di- 
rection is a very strong one, which seems suggestive but has only 
little support from structure formation simulations. In fact, a num- 
ber of studies point at possibly large misalignments and underline 
the complexity of the baryonic physics on galactic scales (van den 
Bosch et al. 2002; Navarro et al. 2004; Bailin et al. 2005; Bailin & 
Steinmetz 2005; Mayer et al. 2008; Kimm et al. 201 1). In our anal- 
ysis, misalignments between the symmetry axis of the galactic disk 
and the angular momentum axis of the host halo could be incorpo- 
rated in choosing a smaller value for the disk thickness parameter a, 
which will play the role of normalising the ellipticity spectra. The 
angular momentum-based alignment model is only able to capture 
the physics of tidal alignment of spiral galaxies. In the case of el- 
liptical galaxies, a model which is linear in the tidal shear is more 
appropriate. 



2.6 Weak gravitational lensing 

The weak lensing convergence k provides a weighted line-of-sight 
measurement of the matter density 6 (for reviews, see Mellier 1999; 
Bartelmann & Schneider 2001; Huterer 2002; Hoekstra & Jain 
2008; Bartelmann 2010) 



k= I d X W K (x)S, 



(12) 



with the weak lensing efficiency W k (x) as the weighting function, 

30 n 4 c XH dz y' - v 

W K (x) = ^—G(x) X , with G( X ) = dx'n(z) — ^-^.(13) 

Jy 



2 xl « 



4r' X' 



n(z) denotes the redshift distribution of the lensed background 
galaxies (with the parameterisation introduced by Smail et al. 
1995), 



n{z) = n 



exp 



dz with 



z _ r (3) 



(14) 



1 

Zo has been chosen to be ^ 0.64 such that the median of the red- 
shift distribution is 0.9, which is anticipated for the EUCLID galaxy 
sample (Amara & Refregier 2007; Amendola et al. 2012). With 
these definitions, one can carry out a Limber-projection (Limber 
1954) of the weak lensing convergence for obtaining the angular 
convergence spectrum C K (l), 
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CM 



Jo x 2 



W 2 ( X )P(k = tlx), 



(15) 



which describes the fluctuation statistics of the convergence field. 
We will always work in the weak lensing regime, k,j <k 1, and 
approximate the reduced shear g = y/(l — k) with the lensing shear 
y, which has the same statistical properties as the weak lensing 
convergence k. 



3 ELLIPTICITY CORRELATIONS 

3.1 Angular momentum induced ellipticity correlations 

The idea behind intrinsic correlations is that neighbouring galaxies 
build up their angular momenta with correlated tidal shears because 
the galaxy separation is typically smaller than the correlation length 
of the tidal shear field. Under the assumption that the galactic disk 
orients itself perpendicular to the angular momentum direction of 
the host halo (for a review on angular momenta of galactic disks, 
see Romanowsky & Fall 2012; Bryan et al. 2012), one perceives 
neighbouring galactic disks under correlated angles of inclination, 
and therefore the apparent shapes are correlated, which is measured 
in terms of ellipticities. We use two ellipticity correlation models 
in this paper, which are both constructed on the idea of correlated 
angular momenta, but which differ in their particular ansatz. The 
first model, proposed by Crittenden et al. (2001) establishes the link 
between the angular momentum direction to the tidal shear field in 
a random process in real space, whereas the second model, which 
is due to Mackey et al. (2002), directly formulates the ellipticity 
field in Fourier space, which makes it easier to quantify ellipticity 
spectra, but whose parameterisation is not as clear as in the first 
case. 



3.1.1 ellipticity correlations 

Correlations of the two ellipticity components e + and e x between 
two points 6] and 2 separated by an angular distance 9 can be 
described using two correlation functions C ++ (0) = (e + (#i)e + (#2)) 
and C xx (#) = (e x (#i) f x(#2))> which are conveniently combined into 
two correlation functions C ± (6), 



C + (8) = C ++ (0) + C xx (9) 
C-(G) = C ++ (0) - C xx (0) 



(16) 
(17) 



using C +x (0) = (e+(#i)e x (02)) = 0. Finally, ellipticity correlation 
functions can be transformed to the spectra C e E (() and C e B {t) of the 
gradient and vorticity modes of the ellipticity field, 



cm 



cm 



n J 9d9[C+(9)J (i9) + C-(9)J 4 (£9)] . 
n J 9d9 [C + ( 



MMm - c^9)j 4 (W)] , 



(18) 



(19) 



by Fourier transform (Kaiser 1992; Schneider et al. 2002; Schnei- 
der & Kilbinger 2007; Fu & Kilbinger 2010). Gravitational lensing 
in the lowest approximation is only able to excite £-modes in the 
ellipticity field. 



3.1.2 configuration space approach 

In this work we use the angular momentum-based ellipticity cor- 
relation model proposed by Crittenden et al. (2001) (referred to 
as the CNPT-model), who trace ellipticity correlations back to 



tidal shear correlations using the conditional probability distribu- 
tion p(L|O tt/3 )dL introduced by Lee & Pen (2001): In this model, 
the distribution p{L\<S> a p)dL is assumed as being Gaussian which is 
then being marginalised over the magnitude of the angular momen- 
tum vector, retaining only its directional dependence. Writing down 
the ellipticity components as a function of the angular momentum 
direction and employing the covariance {L a Lp) as a function of the 
squared tidal shear tensor, as advocated by Lee and Pen, it is pos- 
sible to relate the tidal shear correlations to the spectrum of the 
density field. 

Angular momenta L are described as being coupled to the tidal 
shear by means of a Gaussian random process p(L\<b a p)AL involv- 
ing tidal fields shaping the covariance cov(L) ff/3 of the Gaussian 
distribution (Lee & Pen 2001), 



<L 2 > 



1 +a 



Sa/i - a (U> 2 ) a/ 



(20) 



with the misalignment parameter a, which describes the average 
orientation of the protohalo's inertia to the tidal shear eigensystem. 
a has been measured in numerical simulation to be close to 0.25 
which we will assume in this work. O is the unit normalised trace- 
less tidal shear with the properties tr(O) = and tr(® 2 ) = 1. This 
description is valid on scales where the correlations between inertia 
tensors are negligible. 

The conditional probability density can be used for establish- 
ing a direct relation between ellipticity e and tidal shear <b a p by 
integrating out angular momentum direction and magnitude: 



e(Kfi) = J dLe(L) J L 2 dL p(L\<b a p) 



(21) 



With this relation, one can write down the two correlation func- 
tions (e+(x i)c+(x2)) and (e x (x i)e x (x2)) of the three-dimensional el- 
lipticity field in terms of moments („(r) (see Crittenden et al. 2001) 
of the tidal shear field. Those moments, in turn, are expressed as 
weighted integrals over the CDM-spectrum, where we impose a 
Gaussian smoothing on a scale of 10" M Q /h, which is typical for 
galaxies. 

The correlation function of the 3-dimensional ellipticity field 
can then be projected onto the angular correlation function of the 
ellipticity components by using the configuration-space Limber- 
equation (Limber 1954): 



C ++ (0) = J d Xl W e (xi) J d X 2W € (x2)(e + {xi)e + (x2)) 
C xx (9) = J d Xl W e (Xi) J d^ 2 !V f 0r 2 )(e x (x 1 )e x (x 2 )> 



(22) 



(23) 



with the distance distribution W € (x) = n(z(x))dz/dx resulting for a 
given cosmology from the observed redshift distribution n(z)dz of 
background galaxies (see eqn. 14). The separation distance enter- 
ing the three dimensional correlation functions is completely de- 
termined by the two line-of-sight distances X uX2 an d the angle of 
separation 9. Giahi-Saravani & Schafer (2012) have shown that dis- 
torsions of the intrinsic ellipticity pattern due to the peculiar motion 
of galaxies is very small for multipoles up to t = 1000. 



3.1.3 Fourier-approach 

We extend our analysis by the approach of Mackey et al. (2002) 
(abbreviated as MWK). Similar to Crittenden et al. (2001) they 
also work in the framework of tidal torque theory and relate in- 
trinsic ellipticity to angular momentum assuming that the disk of a 
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galaxy forms perpendicular to its spin axis. However, when com- 
puting the angular momentum the MWK-model entirely neglects 
any correlation between the tidal field and the tensor of inertia. 
They argue that due to the different correlation lengths involved 
(while the correlations for the inertia tensor primarily arise from 
smaller scales the correlations in the tidal field are long-ranged) 
a successive averaging-process is permissible. First, they perform 
an average over inertia tensors then over the tidal field expecting 
the ellipticity correlation arising from long-range correlations of 
the latter. In contrast to this, the CNPT-model takes explicitly the 
correlations of the inertia tensor and the tidal field via the misalign- 
ment parameter a into account. 

One more simplification made by MWK is to drop the depen- 
dence of the observed ellipticity on the z-component of the angular 
momentum of the galaxy, i.e. 



C(l]- Lfj and e x = 2C L X L V . 



(24) 



with a constant C. Accordingly, the intrinsic ellipticity scales 
quadratically with the modulus of the angular momentum, leav- 
ing faster spinning galaxies more flattened. Note that by the basic 
structure of this relation the symmetry properties of the ellipticity 
field are the same ones as in eqn. (1). 

Assuming shear and inertia being statistically independent al- 
lows to work completely in harmonic space, which greatly facil- 
itates the computation. Since we aim at the power spectra of the 
intrinsic ellipticity it is convenient to introduce the parity conserv- 
ing (£-mode) and parity violating (6-mode) part of the intrinsic 
ellipticity field 

E(k) k 2 = [kl - k 2 ) e+(k) + 2k,k y e x (k) 

fi(k)* 2 = -^UW + fe-^f.fk). (25) 



Focusing on modes perpendicular to the line of sight one can derive 
the following dimensionless ellipticity power spectra for the E- and 
6-mode, respectively 



A 2 x (k) = 



C /3 ,\ 4 f°° da . 
225 1 ~~a 



(ak) 



r [ A 2 (/c^l + a--2ap) 
J ^ ,-, . .9 — o_...v7/ 7 Sx(a,n). 



(26) 



(1 + a 2 -2ap) 1 l 2 

Here X e {E, B) and g x is a polynomial given in eqn. (17) of 
Mackey et al. (2002) together with a detailed derivation of the ex- 
pressions given above. As before, we smooth the linear power spec- 
trum on an appropriate length scale with a Gaussian filter function. 

Finally, in order to get the corresponding angular power spec- 
tra we have to again make use of the Fourier-space variant of Lim- 
ber's projection (Limber 1954) 



((2£ + 1) 
Arc 



cm 



X A X Wl( X )A\{llx) 



(27) 



with the weighting function W e C^)d^ = n(z)dz already introduced 
in eqn. (14). 

Mackey et al. (2002) determine the constant C by comput- 
ing the expectation value of the squared angular momentum mod- 
ulus and adjusting C to match the mean-square source ellipticity 
typically observed in galaxy surveys. For our purpose, however, 
it makes more sense to choose C in such a way that the angu- 
lar power spectra obtained with the two different approaches co- 
incide on largest scales. This is justified by the expectation that 
the large-scale power will be least effected by the differences in 
the two approaches under consideration. Hence, we first compute 
the spectra using formulae (18) and (19) for appropriately chosen 




angular scale 8 [arcmin] 



Figure 1. Angular ellipticity correlation functions C++ (9) (green and ma- 
genta lines) and C X x(#) (red and cyan lines), for a smoothing scale of 
M = 10"M G //i, a misalignment parameter a = 0.25 and a disk thickness 
of a = 0.75. The correlation functions were derived using the CNPT- and 
MWK-models with the relative normalisation as discussed in the text. 



misalignment parameter a and galaxy thickness parameter a and 
subsequently determine C, so that the two models yield identical 
predictions for the variance of the intrinsic ellipticity field on large 
scales. 



3.1.4 comparison of the two ellipticity-models 

The ellipticity correlation functions C++{8) and C XX (S) resulting 
from both models are plotted in Fig. 1 as a function of angle of 
separation 6, where the projection was carried out for the EUCLID 
galaxy redshift distribution. The plot suggests correlation lengths 
of ~ 10 arcminutes for the ellipticity field and shows that under 
the normalisations chosen, the correlation functions resulting from 
the CNPT-model achieves 50% higher amplitudes in comparison 
to those predicted by the MWK-model, but otherwise the general 
shape is in very good agreement. 

Computing the spectra C E (£) and C e B (() yields Fig. 2, where 
for comparison the linear and nonlinear spectra C K {() for the weak 
lensing convergence and the EUCLID shape noise levels crl/n are 
plotted. The shape of the ellipticity spectra shows constant ampli- 
tudes up to scales of € ~ 300, where individual ellipticities are 
uncorrelated, and correlations between ellipticities are present on 
smaller angular scales. The spectra exhibit a wide maximum on 
multipoles of I ~ 10 3 before dropping in amplitude, which is 
caused by imposing the mass cutoff. For comparison and motivat- 
ing our analysis we plot predictions for the weak lensing spectrum 
C K (f) for linear and nonlinear CDM spectra, as well as the antici- 
pated shot-noise for EUCLID. Clearly, intrinsic ellipticity correla- 
tions are subdominant compared to weak lensing induced ellipticity 
correlations, but can in amplitude amount to up to 30% of the lens- 
ing signal on multipoles of t - 10 3 before the shape noise makes 
measurements difficult. Comparing the two ellipticity models show 
that, if the normalisation is chosen as explained, the MWK-model 
predicts lower spectra than the model by CNPT by about 50% on 
high multipoles as in the case of the correlation function at small 
separations, but both models predict similar ratios between C e E (t) 
and C B {£) amounting to about a factor of 5 at high multipoles. In- 
terestingly, intrinsic ellipticity correlations would dominate over 
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Figure 2. Ellipticity spectra C E (C) (green and magenta lines) and C B (l) 
(red and cyan lines) derived with the CNPT- and MWK-models, with weak 
convergence power spectrum C K (T), both linear and nonlinear (black dashed 
and solid lines, respectively), and the EUCLID shape noise (r\jh, with 10" X 
<T?/S, n = 0, 1, 2 (blue dashed lines). 



the weak lensing signal if the lensing prediction was derived us- 
ing linear structure growth only. Comparing the spectra C e E {() and 
C K {£) with the shape noise levels of EUCLID clearly demonstrate 
the importance of intrinsic ellipticity correlations in weak lensing 
data. 

The S-mode spectrum C e B (£), which is sourced by intrinsic 
alignments, is smaller by more than one order of magnitude com- 
pared to the £-mode spectrum C e E (£) at high multipoles but might 
dominate over other higher-order lensing effects which are able to 
excite parity-violating modes in the ellipticity field such as source- 
lens clustering (Schneider et al. 2002), multiple lensing along the 
line-of-sight or violations of the Born-approximation (Cooray & 
Hu 2002; Shapiro & Cooray 2006; Schafer et al. 2012). 



3.2 Variance in apertures 

Quantities derived from the weak lensing spectrum C K (f) are 
weighted variances of the convergence inside apertures of varying 
size 6, introduced by Schneider (1996): 



and 



(Mlm = lffM w 2 4 m cao, 



(28) 



(29) 



which measure the scale-dependence of fluctuations' variance in 
the respective field. The weighting functions Wq(x), and VKiOt), x = 
£9, are defined as: 



Figure 3. Standard deviation of the averaged ellipticity and aperture mass 
as a function of aperture size 6 for the CNPT- and MWK-models: tangential 
ellipticity (solid green and magenta lines), radial ellipticity (ej> (solid 
red and cyan lines), aperture mass (M^) (dashed green and magenta lines) 
along with (M^) (dashed red and cyan lines), all for the EUCLID galaxy 
sample in comparison to the same quantities derived from the weak lensing 
convergence (corresponding black lines). 



(M 2 ± xe) 



1 fide Wl(£6)C e E {(), 
^ f £d£ W*(£9)C e B (£), 

m wj(mc e B ({), 



(32) 
(33) 
(34) 



Wi(x) = and W 4 (x) = — , 



(30) 



which are analogous quanitites if the origin of ellipticity correla- 
tions would be purely intrinsic and generated by correlated angular 
momenta. 

In Fig. 3 we show, how the aperture-weighted variances of 
the intrinsic ellipticity field should behave in the EUCLID galaxy 
sample as a function of angular scale in comparison to that of the 
weak lensing field. We consider both tangential and radial shears 
and compare the results between the two different intrinsic ellipitic- 
ity models. 

The aperture-weighted variances derived from intrinsic ellip- 
ticities decrease in magnitude which of course is a generally ex- 
pected behaviour caused by the weighting functions W„((6), and 
exhibit lower amplitudes compared to the lensing ones, as can be 
expected from the relative magnitudes of the spectra. From angu- 
lar scales of 100 arcminutes on intrinsic alignments have dropped 
to zero, which is compatible with them being a small-scale phe- 
nomenon, while weak lensing still has a considerable signal. As 
expected, the differnece in magnitude of the variances sourced by 
£-modes and B-modes is smaller than the difference in spectra on 
small scales because of the averaging in multipole £, and an analo- 
gous argument applies to the predictions by the two models under 
consideration. 



respectively, for the shear variance averaged in an aperture of size 8 
and the aperture mass variance. The weak lensing power spectrum 
is substituted in the previous definitions with C E (£), and C e B (£) in 
order to obtain: 



t&t Wt(£9)C E (£), 



(3D 



4 PARAMETER LIKELIHOOD 

In this section, the dependence of the intrinsic ellipticity spectrum 
on the cosmological parameter set is investigated. This is of partic- 
ular relevance because of their contaminating effect in weak lens- 
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ing data by introducing spurious ellipticity correlations, and be- 
cause they depend on the cosmological model in a very nonlinear 
way, much stronger than e.g. the weak lensing convergence: We 
point out that in our models, the angular momentum L reflects the 
squared tidal shears <9,<9,C>, and the ellipticity field e in turn has 
a very complex dependence on the angular momentum direction, 
which can be approximated to be quadratic for small line of sight- 
components of the angular momentum direction. In the following, 
we keep the parameter a in our ellipticity model constant, because 
it can in principle be determined by analysing morphological data, 
as Crittenden et al. (2001) demonstrated. The misalignment param- 
eter a replaces cr 8 because it fixes the normalisation of the spectra, 
and is completely degenerate with a. 



4.1 Parameter dependences of intrinsic alignments 

When considering the angular spectra describing the ellipticity 
field, the parameter sensitivity of intrinsic alignments can be nicely 
illustrated by considering derivatives of the spectra with respect to 
cosmological parameters, weighted by the inverse noise. For illus- 
tration, we assume that weak lensing-induced ellipticity alignments 
were absent from the data, and that the ellipticity shape noise would 
be that of EUCLID. Therefore the sensitivities, i.e. the derivative 
of the observables with respect to the parameters to be estimated in 
units of the noise, 



1 



dCt 



Vcovx(£) dXf, 



Xe{E,B], 



(35) 



correspond to the contributions -JdF^JdC to the diagonal entries 
of a Fisher-matrix F llv , which describes the parameter dependence 
of the spectra C € E {() and C € B (C) on a cosmological model. The co- 
variances of acquire a cosmic variance error and a Poissonian shape 
measurement error. 



COV X (f) : 



1 



2£ + 1 /. 



sky 



cut) + 



(36) 



with cr E = 0.3, n = 30/arcmin as the number density of galaxies 
per square steradian and the sky fraction / sky = 1 /2. 

We depict these quantities in Fig. 4 for the basic set of cos- 
mological parameters considered here: x fl e {Q. m ,a, h, w), where 
n s has been omitted due to its very weak influence on the spec- 
tra. Clearly, the £-mode and B-mode spectra exhibit an identical 
behaviour on large, cosmic variance dominated scales, where they 
reflect identical dependence on the physical processes of angular 
momentum generation and disk orientation, before differening on 
multipoles £ t 300, where weaker B-modes start being influenced 
by the noise level. The effectively non-existent dependence of the 
ellipticity spectra on the dark energy equation of state parameter 
w is particularly interesting and suggests that intrinsic alignments 
contaminations in weak lensing can be investigated almost inde- 
pendently from the dark energy model assumed. Conversely, the 
dependences on Q„, and a are particularly strong, because a deter- 
mines the amplitude of the spectra in much the same way as erg 
fixes the normalisation of the weak lensing spectra. 



4.2 non-Gaussian likelihoods 

This section is intended to check whether our assumption of Gaus- 
sianity for the parameters likelihood is well-grounded. In order to 
achieve this we compare the Gaussian likelihoods derived by using 
the Fisher-formalism, 



5C E E (l)/Sx^ 

_ac E ji)/ax ii 




100 300 

multipole order I 



1 000 3000 



Figure 4. Sensitivities dC^iO/dx^ (green lines) and BCgiO/dXf, (red lines) 
in units of ycovx(f) as a function of maximum multipole order C, with 
respect to the cosmological parameters il,„ (dashed lines), a (solid lines), 
h (dash-dotted lines) and w (dotted lines). The spectra were computed with 
the CNPT-model. 



£, oc exp 



2cr 2 , 



with cr„ 



(37) 



where the Fisher-matrix is determined from the curvature of the 
logarithmic likelihood, 



dc E (t) i dc E (t) 

dxu cov E (£) dx y 



(38) 



with a direct evaluation of the respective likelihood function, de- 
rived using the relation: 



£ocefsp(-^(^,)/2) 

where the ^-functional is given by 



"max i j 



(39) 



(40) 



and quantifies the goodness-of-fit of a parameter choice x^ 
in the space spanned by the set of cosmological parameters 
jO„,, a, h, w, n s }. It is important to specify that due to the weakness 
of the intrinsic ellipticity correlations we consider only conditional 
errors, i.e. we let one parameter vary at a time, keeping the remain- 
ing fixed and focus only on the stronger £-mode spectrum C e E (C). 
The reduced dimensionality allows us to compute the likelihood 
directly on a grid without having to make use of Monte-Carlo sam- 
pling techniques for evaluating the likelihood. 

Any deviation from a Gaussian shape of the likelihood £. 
is caused by a nonlinear dependence of the spectrum C C E ({) on a 
model parameter x n , which is due to the fact that the x 1 -functional 
deviates from a parabolic shape if the model parameter is varied. 
The likelihood assumes an approximately Gaussian shape if it is 
sufficiently peaked such that a Taylor-expansion of the nonlinear 
parameter dependences is applicable in the region around the fidu- 
cial parameter choice. In our case, non-Gaussian shapes have been 
observed if the summation in eqn. (40) was restricted to the mul- 
tipole range 10 < I < 100 and quickly became Gaussian if the 
summation was carried out to higher multipoles. 
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Figure 5. Conditional likelihoods £.(Sl m ) (solid blue line) and £,(yv) (solid 
red line), along with their Gaussian approximations derived with the Fisher- 
formalism (cr^ = 1/ JFpp, dashed green line), for an observation of elliptic- 
ity spectrum C £ (f) with the EUCLID survey characteristics with all other 
parameters fixed to their fiducial values. The multipole range was set to 
10 < t < 100 and the spectrum C € E (C) entering the likelihood-calculation 
resulted from the CNPT-model. 



The non-Gaussian likelihoods £.{Q.,„) and £.(w) for the mat- 
ter density and the dark energy equation of state can be seen in 
Fig. 5 in comparison to their Gaussian approximation. All likeli- 
hoods are centered on the fiducial model value x^ d and scaled with 
the width cr^ derived with the Fisher formalism. Most notably, the 
likelihoods are more strongly peaked than their Gaussian counter- 
parts, with slight asymmetries of £.{Sl m ) towards large values and 
of X(w) towards small parameter values. It should be noted, that 
the misalignment parameter a (together with the disk thickness a) 
is a linear parameter in our models and its likelihood £.(a) is al- 
ways of Gaussian shape. We conclude that the amount of deviation 
from the ideal shape is not a serious impediment for applying the 
Fisher-formalism for investigating intrinsic alignments, keeping in 
mind that in reality one observes alignments over a much wider 
multipole range such that the likelihoods are closer to Gaussianity. 



5 INTERFERENCE WITH WEAK LENSING 

In this section we show how the intrinsic alignments can affect 
measurements of the convergence spectrum by quantifying the pa- 
rameter biases arising when trying to explain the data consisting 
of both weak lensing-induced ellipticity correlations and intrinsic 
alignments by a model that only accounts for weak lensing and 
neglects intrinsic alignments. We consider the case of EUCLID'S 
weak lensing survey in a non-tomographic setup and give an esti- 
mation of the biases, if intrinsic alignments are not removed from 
data (as proposed by King 2005; Joachimi & Schneider 2009, 2008) 
or not properly modelled (King & Schneider 2003, 2002). We aim 
to supplement previous analysis of intrinsic alignment contamina- 
tions such as Hirata & Seljak (2004), Bridle & King (2007) and 
Kitching et al. (2008) by using a physically motivated and well de- 
scribed alignment model for spiral galaxies with a small number of 
parameters which can be accessed by morphological galaxy sam- 
ples (the disk thickness parameter a) and cosmological simulations 
(the misalignment parameter a). 



5.1 Parameter constraints 

Statistical errors on constraints on the cosmological parameters 
from the projected weak lensing power spectrum C k (l) when ignor- 
ing the intrinsic ellipticity spectrum can be easily obtained by using 
the Fisher matrix formalism, where the Fisher-matrix F^, measures 
the curvature of the logarithmic likelihood In £, in all parameter di- 
rections (Tegmark et al. 1997): 



1 



dC K (0 



1 



3C,(f) 



dxu cow K {() dx v 



(41) 



(42) 



where the covariance is given by: 

2 1/ a- 
coy,(€) = — — - — C,(3 + -£ 
2t + 1 / sky \ n 

In the latter expression EUCLID'S noise a 2 e /n was used as well as 
./sky = 1/2 for the sky fraction, and the parameter space is spanned 
by the cosmological parameters x M 6 {fl,„, cr 8 , h, n s , w}. We will use 
the covariance cov^(^) throughout this chapter and neglect small 
contributions due to intrinsic alignments to the covariance of the 
ellipticity field, as the covariance and therefore the sampling noise 
is dominated by weak lensing. 

The statistical ncr-ellipses obtained as cross-sections through 
the Gaussian-approximated likelihood for all pairs of parameters 
are shown in Fig. 6 together with the systematical errors in param- 
eter estimation if intrinsic alignments are not taken care of. The 
extend to which weak lensing parameter likelihoods are Gaussian 
is investigated in detail by Wolz et al. (2012). 



5.2 Parameter estimation biases 

Ultimately our analysis aims to quantify how biased the parame- 
ter estimation with EUCLID-data will be if intrinsic alignments as 
predicted from angular momentum models were present in the data 
but if we were to interpret the data with a model which does not 
take intrinsic ellipticities into account. We therefore identify a true 
model which includes intrinsic alignments 



C,(f) = C K (l) + C € E (£) + 



(43) 



and a false model, which omits intrinsic alignments and considers 
the ellipticities as random, 



cm = cm + 



(44) 



For both models, one can specify a goodness-of-fit parameter which 
in the case of Gaussian errors is the^ 2 -functional. If the data, which 
follows the model C,(f), is interpreted with the wrong model Cf(f), 
the corresponding^- 2 -functional will exhibit its minimum at a posi- 
tion in parameter space shifted from the true parameter choice, be- 
cause the incomplete model is forced to provide a fit to the data by 
detuning the parameter set away from the fiducial values. The way 
to achieve this goal was proposed by a number of authors (Cabre 
et al. 2007; Amara & Refregier 2008; Taburet et al. 2009; March 
et al. 201 1 ; Schafer & Heisenberg 2012) in different contexts: Tak- 
ing the second-order Taylor expansion of the wrong ^ -functional 
around the best-fit point x, of the true model one retrieves an ex- 
pression involving the vector 8: 

3 1 d~ 

Xj(x f ) = /fix,) + ^ g^xjMS, + - £ ^-q^^,)6,S v , (45) 

ft ' & fi,V y 
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Figure 6. Parameter estimation biases in fi m , o"8, h, n s and w obtained respectively with the CNPT-model (in red) and with the MWK-model (in yellow) for 
EUCLID'S observation of the weak lensing spectrum C r (f), which is contaminated by intrinsic alignments C € E (£) on small scales. The misalignment parameter 
is set to a = 0.25; the disk thickness parameter is set to a = 0.75, the multipole range was 10 < ( < 3000 and as a noise amplitude we considered <x e = 0.3. 
The ellipses give 1 , 2, 3cr statistical uncertainties on the cosmological parameters from the weak lensing spectrum C K (C ) from the same multipole range. 



being 6 = Xf — X t . For weak systematics, this approach has been 
demonstrated to yield very accurate results for the estimation biases 
by comparison with the shift of the likelihood peak evaluated by 
MCMC-techniques (Taburet et al. 2010). 

Now, by extremising the ensamble-averaged (Xf(Xf)) one gets 
the best-fit position Xf. This operation yields a linear system of 
equations: 



6 P = X ( - G '• ) '"' aV 



(46) 



and can be inverted directly leading to the estimation bias 6. Substi- 
tution gives expressions for the quantities G^ v and a p which involve 
derivatives of the spectra, 



Cl " = Z-i covT 



dCMdCM C e (£) d 2 CM 
dxp dx v 

dC K (t) 



dx^dxy 



cm 



(47) 



It is worth to notice that the espression for G,,„ simplifies to F MV 
when the correct model is used, and the bias vector is therefore 
zero. A consequence of this argument is that the inclusion of a 
Gaussian prior FjJ" or would reduce the parameter estimation bias 
due to the transformation G„ v — > G„ v + F^ v l0r , leading to smaller 
values for 6 U in the inversion of the linear system eqn. (46). 

Fig. 6 shows the biases in the estimation of the cosmologi- 
cal parameters induced by considering the intrinsic ellipticities. We 
computed the 1,2, and 3cr ellipses by means of the Fisher matrix, 
as explained in Sect. 5.1, and we considered the full range of mul- 
tipoles going from ( min = 10 to f raax = 3000 well into the noise- 
dominated regime. We calculated the biases for both the CNPT and 
the MWK models of the intrinsic alignments: In what concerns the 
CNPT-model we considered the value of the misalignment param- 
eter found in n-body simulations (Lee & Pen 2000) a = 0.25, and 
a = 0.75 for the thickness of the disk. The biases are shown re- 
spectively in red for the CNPT-model, and in yellow for the MWK- 
model. All parameters apart from the dark energy equation of state 
w and the slope n s are biased significantly by at least lcr, and in 
almost all cases the shift is not along the primary statistical degen- 
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eracy. The differences in the biases between the CNPT- and MWK- 
models reflect the difference in ampiitude they predict, which is in 
our case due to the choice of normalisation. 

The first thing to notice is how the biases depend on the param- 
eters. Evidently Cl m and <x 8 , parameters on which the convergence 
spectrum C K (C) highly depends, seem to be mostly affected, and this 
is in line with the fact that, in a Bayesian analysis, the fit will tweak 
preferably those parameters the model is more strongly dependent 
on, therefore generating larger biases. The estimation biases in Q.„, 
and erg suggest that the presence of C e E (C) is increasing the nor- 
malisation and adds power to the high-f part of the spectrum in 
exactly the way which was expected by choosing higher cr 8 -values 
and lower f2„,-values. For these two cosmological parameters, how- 
ever, excellent priors from the CMB temperature and polarisation 
spectra are available which can be used for further reducing the 
biases. 

Nonetheless our finding of substantial invariance of the dark 
energy parameter w appears to be in contradiction with the results 
obtained by Kirk et al. (2012), who performed a joint analysis on 
galaxy shape and galaxy number density correlations in order to 
set constraints on cosmological parameters, thereby accounting for 
both intrinsic alignments and galaxy bias and found very large bi- 
ases in w. It seems very interesting that the linear alignment model 
used by (Kirk et al. 2012) and our quadratic alignment model lead 
to such diverse conclusions. Due to the complication and difficul- 
ties of physically modelling the intrinsic alignments on the one 
hand, and due to the little information known about the galaxy 
bias on the other, in this work the two contaminants have been 
parametrized as functions of scale and redshift (see also Sect. 2 
of Bridle & King 2007) for further details on intrinsic alignment 
modelling), and their values at different (k, z) points have been con- 
sidered as nuisance parameters along with the cosmological param- 
eters. The constraints on the cosmological parameters have been 
then obtained by marginalizing over the nuisance parameters, and 
the analysis of Joachimi & Bridle (2010) and Kirk et al. (2012) ul- 
timately show that the parameters which are mostly affected by the 
marginalization are the ones entering the dark energy equation of 
state wo and w a . This illustrates the range of predictions from dif- 
ferent alignment models and emphasises the need of a better theo- 
retical understanding of galaxy alignments. 

In fact, this was our motivation for investigating the problem 
of intrinsic alignments by adopting a physical model, although sim- 
plified, which allows to keep the statistics separated from systemat- 
ics, and therefore to quantify the net contribution of the latter to the 
first. On the other hand, the choice of using larger sets of data and 
incorporating the biases by means of a parametrisation certainly in- 
creases the statistics, but necessarily leads to a mixing of statistics 
and systematics, since the biases are disguised as nuisance param- 
eters to be marginalized over. This leads to internally calibrated 
constraints on the cosmological parameters, but makes it difficult 
to estimate to which extent the modelling of the data is actually 
affected by the inclusion of contaminating effects such as intrinsic 
alignments. 

It is worth noting that for linearly evolving Gaussian ran- 
dom fields the quadratic alignment model does not yield a nonzero 
prediction for cross-correlations between intrinsic alignments and 
weak gravitational lensing (Gl-alignments), as this correlation 
would be proportional to the third moment of a Gaussian ran- 
dom field. Only in the case of nonlinearly evolving density fields 
this correlation would be nonzero due to higher-order corrections 
in the weak lensing signal. But conversely, positive detections of 
Gl-alignments on linear scales would be a clear signature of lin- 
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Figure 7. Biases in units of the statistical error cr^ for the cosmological 
parameters Cl m (black line), o"8 (green line), h (red line), n s (magenta line) 
and w (blue line), as a function of the misalignment parameter a of the 
CNPT-model. 

ear alignment models rather than quadratic ones and could help to 
choose the correct family of alignment models. 

5.3 Scaling of the estimation bias 

It is necessary to investigate how the estimation bias scales with the 
normalisation of the intrinsic alignment spectra, as the parameters 
a and a in the CNPT-model and the resulting normalisation C in 
the MWK-model have a large uncertainty. It is worth recalling that 
the misalignment parameter a is measured in ;i-body simulations of 
structure formation, and the galaxy disk thickness a is taken from 
data on galaxy morphologies. 

In Fig. 7 we focus on the CNPT-model, by plotting the biases 
in units of the conditional error a 2 = 1/ F m as a function of the mis- 
alignment parameter a while keeping a fixed at 0.75. An increasing 
value of a means a higher misalignment between shear and inertia 
tensors (Lee & Pen 2000; Crittenden et al. 2001; Schafer 2009), 
and hence a higher correlation between angular momenta, or, dif- 
ferently phrased, less randomness in their directions. This means 
that, within the model, the angular momenta trace the underlying 
tidal shear field in a tighter way. The lower limit a = indicates 
therefore complete randomness and absence of any link of the an- 
gular momenta to the gravitational potential, as described in Lee & 
Pen (2000). It is visible how the biases grow fast for greater values 
of the parameter a, simply meaning that stronger intrinsic correla- 
tions would represent a stronger contamination to the covergence 
spectrum, but that for very large values of a the dependence of the 
parameter estimation bias with a is saturated and it evolves weaker 
with increasing a, and even decreases in the case of the parame- 
ter n s . Given the fact that the estimation biases 6 M / 'cr^ change by 
almost two decades as a is varied, it is vitally important to deter- 
mine a beforehand, either from independent observations or from 
simulations of structure formation. 

5.4 Observations of intrinsic alignments 

Another approach to gain infomation about the intrinsic alignments 
relies on the subtraction of the weak lensing signal from the over- 
all signal that is measured. This is, roughly speaking, the inverse 
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of what is usually done with nulling techniques (for details see 
Joachimi & Schneider 2008, 2009, 2010), whose aim is to clear up 
the weak lensing signal from intrinsic alignments. The idea behind 
this is that if we know with high accuracy the model describing 
the cosmology, we can then predict how, according to this model, 
the weak lensing spectrum must be. By subtracting the latter from 
the measured spectrum, the remaining part is ascribable to intrinsic 
ellipticities. More precisely, if the uncertainty on the covergence 
power spectrum is small enough to still allow for the extraction of 
the intrinsic ellipticity signal, then it means that we will be able, 
in future surveys, to notice this weak signal in the presence of the 
much stronger weak lensing spectrum. 

It is worth to remark, at this point, that the intrinsic align- 
ments we consider in this work, also known as II (intrinsic- 
intrinsic) alignments, are not the only contaminant to the weak 
lensing spectrum, usually referred to as the GG (gravitational shear- 
gravitational shear) signal. For instance, another source of contam- 
ination are the GI (gravitational shear-intrinsic ellipticity) correla- 
tions, which might occur when the alignment produced by a dark 
matter halo on a closeby galaxy correlates with the shear signal that 
the same halo induces on a background galaxy. This effect, first 
suggested by Hirata & Seljak (2004), is difficult to remove, and is 
not considered in our treatment. Likewise, we do not consider com- 
plications arising from the statistical uncertainty in estimating the 
combined spectrum CM) + C e E {£). 

The problem is therefore now to understand whether the un- 
certainty at which the lensing spectrum can be predicted for a given 
cosmology is enough to attribute a high-f excess in the ellipticity 
spectrum to intrinsic alignments. In order to quantify the uncer- 
tainty on the covergence spectrum AC K (€), we consider a multi- 
variate Gaussian likelihood for the cosmological parameters, draw 
from this distribution simultaneously a sample of five parameters, 
and compute for those the weak lensing spectrum. This gives us a 
bundle of spectra around the fiducial spectrum C k .({\ x ff), and allows 
us to define the uncertainty as the standard deviation: 

ACA ° 2 = ~ E [ C M&> - <C*(^ d )>f , (48) 

at each multipole £ where the index i runs over the samples of 
parameter sets drawn from the multivariate Gaussian likelihood. In 
short, this sampling of a parameter set and measuring the variance 
of the resulting spectra is a method of propagating the statistical 
parameter uncertainty described by the likelihood to an error-tube 
around C K (C) reflecting the prediction uncertainty in the spectrum. 
This allows now to compare the magnitude of the intrinsic elliptic- 
ity spectrum C E {1) to this uncertainty and to quantify the signifi- 
cance. 

Fig. 8 shows the uncertainties n x AC K {€) with » = 1 ■ ■ ■ 5 ob- 
tained by using this technique. The likelihood from which samples 
on the wCDM-parameter set including Q m , erg, h, n s and w were 
drawn is the one describing the knowledge on the cosmological 
parameters if EUCLID'S measurement of baryon acoustic oscilla- 
tions, EUCLID's weak lensing data (both in a 10-bin tomographic 
measurement) and PLANCK's CMB data were present, with a 
theoretical prior on spatial flatness. The PLANCK-likelihood was 
marginalised over the baryon density £l b . The error tube AC,(f) 
can be clearly separated into two multipole ranges, the first region 
where the linear CDM-spectrum is dominating in the generation of 
the weak lensing spectrum, and the second region where the non- 
linear enhancement of P(k) is important and where the error tube is 
much wider. 
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Figure 8. Uncertainties AC K (C) (shaded area) in the prediction of the non- 
linear weak lensing spectra C K (-6) from drawing samples for Sl m , erg, h, n s 
and w from a Gaussian parameter likelihood, for which we use a prior on 
the wCDM-model combining baryon acoustic oscillations, lensing and the 
CMB (EUCLID 10-bin BAO spectra, EUCLID 10-bin weak lensing spectra 
and PLANCK CMB temperature and polarisation spectra). The uncertainty 
is compared to the ellipticity spectra C e E {t) (green line) and C e B (l) (red line) 
determined with the CNPT-model. 

With such small uncertainties and with the choice of a cos- 
mological model with low complexity the predictive uncertainty 
on C K (£) is much smaller than the amplitude of the intrinsic align- 
ments from multipoles of I =± 30 on, AC K (£) <k C e E (C). We verified 
that the PLANCK CMB likelihood alone would not be sufficient 
for extracting the intrinsic ellipticity spectrum. Likewise, a more 
complex model with a larger number of parameters would have 
much larger uncertainties, AC„({) ^ C E (£) for most of the multi- 
pole range. 



6 SUMMARY 

Subject of this paper are the statistical properties of intrinsic, angu- 
lar momentum induced ellipticity alignments, and their dependence 
on the cosmological parameter set, in comparison to ellipticity cor- 
relations induced by weak gravitational lensing. We carry out our 
computations with the EUCLID ellipticity data sample in mind, and 
use the projected EUCLID galaxy redshift distribution and shape 
noise for making forecasts. 

(i) We base our predictions for the spectra C e E (E) and C B {£) de- 
scribing fluctuations in the ellipticity field on physical models for 
angular momentum correlations in the large-scale structure (Crit- 
tenden et al. 2001; Mackey et al. 2002). The two models under con- 
sideration link the angular momentum field to the tidal shear field, 
and model the ellipticity of a galaxy by assuming that the galactic 
disk is formed perpendicular to the host halo's angular momentum 
direction. The two models differ in describing these physical pro- 
cesses in configuration space versus Fourier space, and use different 
normalisations. For comparability, we have normalised the MWK- 
model such that it displays the same amplitudes as the CNPT-model 
on large angular scales. The CNPT-model in turn uses 3 parameters, 
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which are the mass-scale of the galaxies, imposed by an Gaussian 
filter acting on the CDM-spectrum P(k), a misalignment parameter 
a, which is determined to have the numerical value a = 0.25 in 
numerical simulations (Lcc & Pen 2001) and finally the disk thick- 
ness parameter a, which has been measured to be a a 0.75 in the 
APM-galaxy sample (Crittenden et al. 2001). 

(ii) Computing the ellipticity spectra C E {£) and C e B (C) for the 
galaxy sample of EUCLID from both models yields spectra which 
are constant on large angular scales and drop off exponentially on 
small scales, while the E-mode spectrum is larger by about an or- 
der of magnitude compared to the S-mode spectrum on multipoles 
of £ 1000. By themselves, the spectrum C E (£) would be signif- 
icantly larger than the linear weak lensing convergence spectrum 
C K (£), which is comparable in amplitude to the spectrum C e B (£) on 
these multipoles. Nonlinear structure formation, however, increases 
the variance of the cosmic density field strongly, such that intrin- 
sic ellipticities contribute only ~ 20% to the total variance of the 
ellipticity field at £ = 1000. Aperture weighted variances give a 
similar impression: The averaged shear and the aperture mass of 
the nonlinear weak lensing convergence dominate the variance on 
all relevant scales, and intrinsic alignments are smaller by at least a 
factor of two in this observable. 

(iii) Investigating the dependence of the ellipticity spectra on 
cosmological parameters gives a result very different compared to 
other cosmological probes. Due to the dependence of the angular 
momentum field on the angular momentum direction and not the 
magnitude, cr 8 is entirely replaced by the misalignment parame- 
ter a. n s , Q.,,, and h determine the CDM spectrum P(k) (the latter 
two by fixing the shape parameter T) and Q.,„ is of course appear- 
ing in the conversion between comoving distance and redshift at 
the stage of applying the Limber equation. Computing the deriva- 
tives dCpify/dXf, and dC B (£)/dx fl of the spectra with respect to the 
cosmological parameters and expressed in units of their covariance 
suggests that the parameters a and fl,„ are the ones most important 
for intrinsic alignments, with only minor dependences on the dark 
energy equation of state w and the Hubble-parameter h. At the same 
time fi„, and tr 8 are the ones best constrained by lensing, so that it 
is suggestive to expect the largest estimation biases in those two 
parameters, if the intrinsic alignments are not properly removed or 
modelled. 

(iv) In the next step we quantified the likelihood £,(Q. m ,a, w) 
of the intrinsic alignment spectrum C e E (£) if lensing was not 
present. One could expect this likelihood to have non-Gaussian 
contributions because of the nonlinearities present in an angular 
momentum-based alignment model: firstly, the angular momentum 
depends on the quadratic tidal shear and the ellipticity depends on 
the squared angular momentum direction. For a resticted multipole 
fmax ~ 100 range one can see clear deviations from Gaussianity, 
which quickly vanish if the multipole range is extended. The mis- 
alignment parameter a, which describes the normalisation of the 
spectra, is described by a Gaussian likelihood and enters as a pref- 
actor. 

(v) We compute estimation biases on the wCDM parameter set 
if intrinsic alignments are not removed from the ellipticity spec- 
trum, i.e. if the data is in reality described by C K (£) + C e E (£) and 
wrongly fitted by C K (£) only. The strongest biases are present in 
Q.,„, which is measured to low and cr 8 , which is estimated too high, 
both at the level of ~ 2<x, making the estimation biases significant. 
Interestingly, the dark energy equation of state is almost unbiased, 
indicating that dark energy investigations are not directly affected 
by intrinsic alignments. Changing the magnitude of the intrinsic 
alignment contamination by increasing the misalignment parame- 



ter a shows an monotonic increase of the estimation biases for all 
parameters except n s where the estimation bias saturates at a m 0.5 
and drops for higher amplitudes. Clearly, these results demand a 
good external prior on a, either from independent measurements or 
from numerical simulations. 

(vi) Finally we investigate if the weak lensing convergence 
spectrum C K {£) can be predicted precisely enough such that a de- 
viation can be attributed to a contribution C E (C). For this purpose, 
we develop a technique for propagating the uncertainty in the set of 
cosmological parameters to the variance AC,(£) 2 around C K {£) for 
the fiducial cosmology. Comparing this uncertainty with the am- 
plitudes C E (f) suggests that it should be measurable at high multi- 
poles. In this process we used a Gaussian likelihood £, on a standard 
w CDM-cosmology reflecting the knowledge on the cosmological 
parameters from EUCLID'S BAO- and weak lensing spectra and 
from the temperature and polarisation spectra of the cosmic mi- 
crowave background measured by PLANCK. 

We plan to extend our research to the intrinsic alignment con- 
tamination of tomographic weak lensing data, and to include cross- 
correlations between the weak lensing shear and the intrinsic el- 
lipticity field, the so-called Gl-alignments, which we aim to de- 
rive from angular momentum-based alignment models and which 
enter the ellipticity spectra at higher order. These Gl-alignments 
are challenging to describe as they introduce ellipticity correlations 
across tomography bins. Additionally, we aim to include a linear 
alignment model for elliptical galaxies and to work with a proper 
morphological mix of ellipticities, working towards a more com- 
plete physical description of alignments in tomographic weak lens- 
ing data. 
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